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WEIGHING THE NON-TRANSITING HOT JUPITER r BOO B 

F. Rodler 1 , M. Lopez-Morales 1-2 , and I. Ribas 1 
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ABSTRACT 

We report the detection of the orbital velocity of non-transiting hot Jupiter r Boo b. By employing 
high- resolution ground-based spectroscopy around 2.3 /im during one half night, we are able to detect 
carbon monoxide absorption lines produced in the planet atmosphere, which shift significantly in 
wavelength during the course of the observations due to the orbital motion of the planet. This detection 
of the planetary signal results in the determination of the orbital inclination being i = 47tg degrees 
and furthermore allow us to solve for the exact planetary mass being m p = 5.6 ± 0.7 Mj up . This 
clearly confirms the planetary nature of the non-transiting companion to r Boo. 
Subject headings: planets and satellites: atmospheres — planets and satellites: individual (r Boo b) 
— techniques: radial velocities — techniques: spectroscopic 
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1. INTRODUCTION 

The radial velocity (RV) technique is the most suc- 
cessful exoplanet detection method, with more than 600 
planet^ detected to date. This technique monitors the 
variation of the RV of the star due to the gravitational 
pull of an unseen companion. Owing that the RV of 
the planet is not measured directly, neither the orbital 
inclination i of the planet nor its exact mass can be de- 
termined. Instead, only an estimation of the minimum 
mass of the planet can be derived by applying Kepler's 
Laws, which give m Pjm i n = m p sini, where m p is the un- 
known planetary mass. For transiting planets - i.e. for 
planets which periodically occult their host stars - the 
orbital inclination i and therefore the exact planetary 
mass can be measured. However, the vast majority of 
extrasolar planets found to date are non-transiting and 
therefore without an exactly determined mass. Conse- 
quently, many of these non-transiting planets, especially 
the ones with largest minimum masses, can only be la- 
belled as planetary candidates, owing to the case that for 
very low orbital inclination they might turn out to be 
brown dwarfs. 

One strategy to determine the exact mass of non- 
transiting planets is to directly measure the planetary 
RV signal via high-resolution (i.e. R = A/AA > 40, 000; 
A denotes the wavelength) spectroscopy with very large 
telescopes. Key to this method is to observe a large num- 
ber of spectral features coming from the planet and to 
observe them at different orbital phases so the traveling 
faint planetary signal can be disentangled from the dom- 
inating stellar one. In the past, several high-resolution 
spectroscopy campaigns at the optical with the goal to 
detect starlight reflected from hot Jupitcrs (i.e. massive 
planets that are a few stellar radii away from their host 
stars) and to measure their exact masses were carried 
out. Although all of them resulted in non-detections 



of reflected light, these campaigns confirmed the low 
reflectivity of hot Jupite rs at visual wavelength s (e.g . 
Charbonneau et ail 119991: iCollier Cameron et al.l l2002t 
Leigh et al.l 120031: I Rodler et all 120011 Towards near- 
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infrared (NIR) wavelengths, the planet-to-star flux ratios 
drastically incre ase due to the strong the r mal emission 
of hot Jupit e rs. Wiedemann et al.l (120011): IBarnes et al.l 
((20071 l200i I2010D and ICubillos et alffeoilh observed 
hot Jupiters by means of high-resolution spectroscopy at 
near-infrared wavelengths, but were not able to detect 
any molecules in their atmospheres. 

The molecule carbon monoxide (CO) is one key to 
detect the radial velocity of an exoplanet, since it ex- 
hibits a dense forest of deep absorption lines in a spec- 
tral band around 2.3 /an. Models predict CO to be 
one of the most abundant molecules in hot gas giant 
exopl anets (jCooper fc Showmanll2006t iSharp fc Burrows! 
|2007|) . To test those models, some papers have already 
been published reporting the detection or suggesting the 
presence of C O in the atmosphere of some transiting 
hot Jupiters (IGrillmair et al.l l2008t ISwain et al.l 120091: 
iDesert et ail 120091: iTinetti et alj|2010l) . All these mea- 
surements were obtained with very low-resolution (R < 
40) HST NICMOS NIR spectra or broad Spitzer pho- 
tometry beween 3.6 and 24 /im, both via transmission 
(primary transit) or emission (secondary eclipse) obser- 
vations. ISnellen et all d20loh reported the detection of 
CO via the analysis of the transmission spectrum of the 
transiting planet HD 209458b using high-resolution spec- 
troscopy between 2.30 and 2.33 /im. This allowed these 
authors to directly measure the RV of a transiting hot 
Jupiter for the very first time. We note that all those 
CO detections were done on transiting planets, but nei- 
ther atmospheric chemicals nor RV shifts have yet been 
detected for non-transiting exoplanets. 

Here we present the results of our attempt to detect 
the motion of the non-transiting planet r Boo b via ob- 
servations of its atmospheric CO spectrum in the NIR, 
by using a technique similar to ISnellen et al.l ([2010). The 
main goal of our study has been the measurement of the 
orbital inclination and the absolute mass of the planet, 
but as a by-product we can also confirm the presence of 
CO. 

2. THE PLANETARY SYSTEM OF r BOO 
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TABLE 1 

Parameters of the star t Boo and its planetary 
comp anion. Abbreviations for the references are: B06 = 

IBUTLER ET AL.I (120061) AND REFERE NCES THERE I N, 

B97 = I Baliunas et al.1 (119971) , G98 = IGonzalezI lfl998l) 
HOP = IHenry et al I 1120001) , VF05 = IValenti fc FischerI 
120051) , VV09 = Ivan Belle ik. von BraunI 11200911 . 



Parameter 


Value 


Error 


Ref. 


Star: 








Spectral type 


F7 IV-V 




B97 


K (mag) 


3.36 


0.05 


VV09 


m* (M Q ) 


1.33 


0.11 


VF05 


R* (R©) 


1.31 


0.06 


G98 


T eS (K) 


6360 


80 


G98 


Frot (d) 


3.31 




B97 


v sini (km s —1 ) 


14.9 


0.5 


H00 


Age (Gyr) 


1.3 


0.4 


VF05 


Planet: 








m p sini (M Jup ) 


4.1 


0.34 


B06 


a (AU) 


0.0481 


0.0028 


B06 


e 


0.023 


0.015 


B06 


K* (km s" 1 ) 


0.4611 


0.0076 


B06 


Orbital period P or b (d) 


3.312458 


0.00002 


new 


T^o (JD) 


2 454 267.497 


0.0122 


new 



Tau Boo b is a massive hot Jupiter, orbiting its host 
star every 3.31 days, and is estimated to have an at- 
mospheric ternjjera ; ture_ab^ve_J_800 K fol lowing the for- 
malism in lLopez-Morales fc Seagerl (|2007f ) assuming zero 
albedo and no energy redistribution from day-to-night 
side. Table [T] summarizes the parameters of the planet 
and its very bright host star. We updated the orbital 
ephemerides by an RV-analys is of high-reso l ution UVES 
spectra of r Boo published in lRodler et al.1 (|2010f ) 

The maximum possible RV semi-amplitude of the 
companion t Boo b is calculated to be -Kp iinax = 
157.3 ± 4.4 km s" 1 via 

^ -Torb ' 

where G is the gravitational constant, the stellar 
mass, and P rb the orbital period. Due to the absence 
of transits in high-precision photometry (Henry et al. 
2000), orbital inclinations larger than i = 83° can be 
excluded. Baliunas et al. (1997) found that the star 
r Boo rotates rapidly with a period commensurate with 
the orbital period of th e pl anet, suggesting tid al lock- 
ing. iCatala et all (|2007t ) and lDonati et all (|2008l) carried 
out spectropolamctric observations of r Boo and mea- 
sured differential rotation in the star, confirming that its 
rotation is synchronized with the orbital motion of the 
planet and s uggest ing an orbital inclination around 40°. 
iLeigh et all (|2003l) and Rodler et al. (2010) attempted 
to measure starlight reflected from r Boo b, and found 
candidate features of marginal significance, indicating or- 
bital inclinations of i — 37° and 41°, respectively. Un- 
der the assumption of a tidal lock and that the stellar 
equator and the orbital plane are co-aligned, Rodler et 
al. (2010) predicted an orbital inclination of i m 46° and 
a planetary mass of m p « 5.6 Mj up . 

3. OBSERVATIONS AND DATA REDUCTION 



We observed r B oo during 5 hours on June 10, 2011, 
by using CRIRES (|Kaeufl et al.ll2004D . mounted on the 
VLT/UT1 on Cerro Paranal, Chile. We collected a total 
of 116 high- resolution spectra of our target plus 2 spec- 
tra of the A-type star HD 116160, which does not show 
any intrinsic absorption lines in the observed wavelength 
regime and therefore could be used for the calibration 
of our telluric model. The date of the observations was 
selected in such a way that the observations were carried 
out at orbital phases between </> = 0.55 and 0.61 (</> = 
corresponds to the mid-transit position if i w 90°), when 
the hot day side of the planet facing the star was largely 
visible and consequently appeared bright in the NIR. In 
addition, based on the ephemerides in Table 1 and the 
value of K+, we expected a variation of the RV of the 
companion up to 40 km s _1 during the course of the 
observations. 

The observations were carried out with the CRIRES 
standard setting of order 24 and a reference wavelength 
of A = 2.3252 /j,m. We used the 0.2" slit to achieve a 
spectral resolution of R = 100,000. We made use of 
the AO-system of CRIRES fo minimize slit losses. Ob- 
serving conditions were good, with a seeing between 0.6" 
and 1.8". The observations were taken at two different 
nodding positions along the slit with the intention to re- 
move the faint OH-lines in the sky background. We chose 
the integration times in such a way that the peak count 
rates per exposure were low (~ 4000 counts) to avoid 
non-linear sensitivity in the four detectors. 

The data were reduced with IRAlfl. All four detec- 
tors showed a pattern aligned with the reading direction. 
This pattern consisted of alternating rows or columns of 
larger and smaller intensities than the mean value and 
is called odd-even effect. We applied a correction for 
the odd-even effec10, which was most present in the data 
recorded with detectors 1 and 4, for which the reading 
direction was perpendicular to the spectral dispersion. 
Since after this step strong noise features coming from 
those two detectors were still present in the data, we 
discarded those two detectors from any further analy- 
sis. The frames were taken in an A-A-B-B-B-B-A-A- 
sequence, where A and B denote the different nodding 
positions along the slit. For all frames as well as for each 
of each two remaining detectors, we combined the two 
consecutive frames, which had been taken at the same 
nodding position. This reduced the number of spectra to 
58. Nodded images were then subtracted to remove sky 
background and dark current. White-light spectra ob- 
tained with the same instrumental configuration in the 
afternoon before and after the observations were used 
to flat-field the data. Using the apall task, we first 
identified and optimally centered the orders in the two 
individual nodding frames and traced these orders by 
adopting a second-order Legendre polynomial along the 
dispersion axis. We then extracted the one-dimensional 
spectra in detectors 2 and 3 and calculated a second- 
order polynomial wavelength solution by adopting as a 
reference system the telluric lines present in the spec- 
tra. The vacuum wavelengths of the telluric lines were 
identified using the line-by-line radiative transfer model 

4 IRAF Project Home Page: |http: / /iraf.noao.edu/| 

5 An algorithm to correct this effect is provided at 
http: / / www.cso.org/sci / facilities / paranal / instruments / crires / tools / 
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(LBLRTM ) routine, which is based on the FASCODE 
algorithm (jClough et al.l 11992ft . The detectors 2 and 3 
covered the wavelength regions 2.303 to 2.317 /jm and 
2.319 to 2.332 /xm, respectively, with a pixel size corre- 
sponding to w 1.6 km s . We furthermore identified 
and removed obviously bad pixel in each spectrum com- 
ing from defects of the detector. 

4. DATA ANALYSIS 

4.1. Overview 

The observed target spectra were heavily contaminated 
by telluric lines, and to much less extent, by the rotation- 
ally broadened stellar absorption lines of the host star. 
A crucial step of the data analysis was the removal of 
the telluric and stellar spectra, and the search for the 
planetary spectrum in the residuals. In the following, we 
provide a description of the different data analysis steps. 

4.2. Step 1: Determination of the instrumental profile 

To calculate the instrumental profile (IP) of the spec- 
trograph for each target observation, we made use of the 
telluric contamination, which dominated our data, and 
of the telluric model spectra (see Step 2). In contrast to 
the large number of telluric lines, the stellar spectrum of 
r Boo exhibits very few absorption features in the ob- 
served wavelength regime, which appeared rotationally 
broadened due to the high wsini (Table 1). Those re- 
gions of stellar absorption features were excluded from 
the determination of the instrumental profile. We nor- 
malized the spectra in such a way that the flux in the 
telluric/pseudo-stellar continuum was at one, and stored 
the photon noise error per spectral pixel. We then de- 
termined the IP as the sum of se ven Gaussian profile s 
in a similar way as described in IValenti et al.l (|l995): 
Around a central Gaussian we grouped three satellite 
Gaussians on each side of it, which allowed us to ac- 
count for asymmetries in the IP. Free parameters were 
the width of the central Gaussian, plus the amplitudes 
of the six satellite Gaussians, while the positions and the 
widths of those satellites were fixed and set a priori in a 
way that their half-widths overlapped. By employing a 
Brent algorithm, we finally determined the free parame- 
ters and convolved the telluric model spectrum with the 
IP. We found that for all spectra, the shape of the IP 
remained constant. However, the width of the IP was 
slightly different at each nodding position and was also 
changing during the course of the night. 

4.3. Step 2: Telluric model 

For the calculation of the atmospheric transmission 
spectrum, we used the LBLRTM code, which is avail- 
able as a Fortran s ource cod As molec ular database we 
adopted HITRAN (jRothman et al.ll2005D . which contains 
the 42 most prominent molecules and isotopes present in 
the atmosph ere of the Earth. Fo llowing the approach 
presented by iSeifahrt et"aT1 (|2010f ). we created a high- 
resolution theoretical vacuum-wavelength telluric spec- 
trum for each observed spectrum by accounting for the 
air mass of the star as well as the weather conditions 
(water vapour density column, temperature and pres- 
sure profiles) during the observations. We retrieved the 

6 Source code and manuals are av ailable under 
http : //rtweb . aer . com/lblrtm_description . html 



weather information from the Global Data Assimilation 
System (GDAS). GDAS models are available in 3 hour 
intervals for any location around the glob^H- 

We first calibrated the line depths of the telluric lines 
by using the A-star observations, which had been taken 
at the beginning of the night. For five water and methane 
lines, the line depths were systematically underestimated 
in the HITRAN model; for these lines we consequently 
updated the abundances in our telluric model. In ad- 
dition to that, the amount of water vapor in the at- 
mosphere constantly changed during the observations, 
which resulted in variations of the line depths of the wa- 
ter lines in the observed spectra. Consequently, for each 
observed target spectrum, we created a set of two tel- 
luric model spectra, one for water and one for the other 
molecules, most notably CH4. We then globally scaled 
the line depths for those two models for each target ob- 
servation. Finally, we multiplied the two different high 
resolution telluric spectra, convolved it with the IP to 
create the telluric model. 

4.4. Step 3: Stellar model 

The few stellar absorption lines remained in the 
data after the removal of the telluric lines. We at- 
tempted to subtract the stellar spectra from the data 
using models, bu t found that existing models (MARCS, 
Gustafsson etaTl |2008; PHOENIX, lHauschildt eTabl 
1991 " do not reproduce the observed stellar spectrum 
of r Boo. The best alternative was to create a high 
signal-to-noise ratio template by co-adding all the ob- 
served telluric-free spectra. Due to the expected large 
orbital motion of the planet, contributions to this tem- 
plate spectrum coming from the planet, which is esti- 
mated to be a thousand times fainter than the star, were 
smeared out to a large extent. Before combining the in- 
dividual stellar spectra to one template spect rum, we ap- 
plied a Savitzky-Golay smoothing algorithm (jPress et al.1 
1992) to each of them, which ensured that the broad stel- 
lar lines remained in the spectra, while the sharp and 
unseen planetary lines were smoothed out. In the final 
step, we subtracted the stellar template spectrum from 
all telluric-free target spectra. 

4.5. Step 4 : Searching for the CO features 

The mean RV of the host star r Boo with re- 
spect to the barycenter of the Solar system is —15.8 ± 
0.5kms _1 (blue-shifted: iBarbier-Brossat & Figonl2000f) . 
The barycentric RV of the Earth with respect to the sky 
position of the target at the time of observations was be- 
tween -22.8 km s^ 1 and -23.4 km s _1 . This means that 
during the times of our measurements, the stellar spec- 
trum of t Boo was observed with a red-shift of about 
8 km s _1 . That red-shift was subtracted from the resul- 
tant RVs of the planetary signal, as described below. 

After the removal of the stellar and telluric absorp- 
tion lines, we corrected the residual spectra for trends 
originating from intra-pixel variations. We then cross- 
correlated the telluric-free and stellar-free target residual 
spectra with a CO-model spectrum (F ig. 1), which was 
calculated with the PHOENIX code (lHauschildt et al.l 
[19971 IHelling et al.ll2lM IWitte et al.ll2011h for a brown 

7 GDAS webpage: http://ready.arl.noaa.gov/READYamet.php 
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dwarf having a temperature of T c g = 1800 K and So- 
lar metallicity. To this end, we shifted the CO-model 
spectrum for each residual spectrum according to the in- 
stantaneous radial velocity V p ((j>,i) with respect to the 
star, which depends on both the orbital phase <f>, which 
was a priori known, and the unknown orbital inclination 
i as: 

V p = -Kp im axSin27r0sini = K p sin 2tt<P . (2) 

For all residual spectra, we then calculated the cor- 
relation value for different RV semi-amplitudes K p of 
the planetary signal and finally determined the cross- 
correlation function (CCF) with respect to K p . 

4.6. Step 5: Confidence level 

The confidence level of the strongest peak of the CCF 
was determine d by employing a boot strap randomisation 
method (e.g. iKuerster et al.lfl997f l: We assigned ran- 
dom values of the orbital phases to the observed spec- 
tra, thereby creating N different data sets. Any signal 
present in the original data was then scrambled in these 
artificial data sets. For all these randomized data sets, 
we again evaluated the model for the free parameter K p 
and located the best fit with its specific CCF-peak value. 

The confidence level of the CCF peak was estimated 
to be ~ 1 — m/N, where m is the number of the best fit 
models having a CCF-peak value larger or equal than the 
maximum of the CCF found in the original, unscrambled 
data sets. Notice that we consider a detection to have a 
confidence of > 0.9987, which translates to > 3a. 

5. RESULTS AND DISCUSSION 

As illustrated in Figures 2 and 3, the results of our 
cross-correlation analysis using a CO-model spectrum 
reveal Doppler velocity shifts consistent with a maxi- 
mum RV semi-amplitude for the planet of K p = 115 ± 
11 km s _1 . The lc-errorbar contains both the actual 
error of the velocity meas urement, which was determined 
via bootstrap resampling (|Barrow et alj|1984 ). as well as 
the shift in velocity due to the uncertainties in the orbital 
solution. In addition to that, our bootstrap randomisa- 
tion analysis revealed that this signal is significant at the 
3.4cr confidence level, with 2 false positives in 3000 trials. 

The measured RV semi-amplitude corresponds to an 
orbital inclination for the system of i = 47lg degrees. 
Adopting the measured orbital inclination, the value of 
the stellar mass in Tablel and Kepler Laws of planetary 
motion, we finally are able to derive an absolute mass for 
the planet of m p = 5.6±0.7 Mj up . This value clearly con- 
firms the planet hypothesis for the hot Jupiter r Boo b. 

Given the precision of our data, possible residuals of 
the planetary spectrum in the stellar template, the lack of 
detailed atmospheric temperature-pressure profile mod- 
els for this planet, and uncertainties in the possible level 
of masking between CO and CH4 (b oth molecules are 
expec ted to co-exist in r Boo, e.g. iBurrows fc Shard 
1999), we cannot provide a reliable estimation of the CO 
abundance in the atmosphere of r Boo b. However, we 
can place a lower limit to the planet-to-star flux ratio in 
the observed wavelength range based on the atmospheric 
model we adopted. To this end, we created data sets 
which were based on our residual spectra. For each of 
these data sets, we injected an artificial planetary CO 
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Fig. 1. — Bottom: Residual spectrum after subtraction of the 
stellar and telluric absorption lines from the data. The upper 
spectrum depicts the model of the dense forest of carbon-monoxide 
(CO) absorption lines. For clarity, the line depths in the CO model 
are shown for a planet-to-star flux ratio of 2 X 10 -3 , and the spec- 
trum is shifted up by 3 units. 
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Fig. 2. — CO absorption in the planet atmosphere of r Boo b. 
The individual cross-correlation functions (CCFs) of the 58 resid- 
ual spectra with the CO model spectrum in the rest-frame of the 
star are shown. During the course of the observations, the or- 
bital motion of the planet produces a RV-shift starting at about 
-35 km s" 1 and ending at 65 km s — 1 , respectively for orbital phases 
of 0.55 and 0.61. The linear grey-scales indicate the strength of 
the cross-correlation signal (bright means absorption, dark means 
emission). 

signal having a RV semi-amplitude of K p = 115 km s _1 
with selected planet-to-star flux ratios from 10 -2 down 
to 10 -4 . To avoid an overlapping of the CCF peaks with 
the real signal, we red-shifted the spectrum of the in- 
jected planet by 50 km s^ 1 . For each data set, we then 
determined that peak of the CCF apart from the one 
of real signal and determined its confidence level. We 
find that with our data set and our CO-model, we are 
sensitive to detect a planet with planet-to-star flux ratio 
larger than m 7 x 10 -4 at the 3c confidence level. 

We emphasize that this measurement of the exact mass 
of t Boo b represents the first successful determination 
of the mass of a non-transiting planet by means of high- 
resolution spectroscopy. This technique has therefore the 
potential of providing direct masses and estimation of the 
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Fig. 3. — The CO signal co-added for all spectra. The cross- 
correlation functions (CCF) from all 58 residual spectra were com- 
bined assuming different RV semi-amplitudes of r Boo b. The peak 
of the CCF occurs at a RV semi-amplitude of 115 ± 11 km s , 
which corresponds to an orbital inclination of i = 47^1 degrees. 
The peak of the CCF is significant at the 3.4c confidence level. 



atmospheric composition of non-transiting exoplanets in 
the near-future, in particular as new facilities like E-ELT 
become available. 



Based on observations made with ESO Telescopes at 
the Paranal Observatory under programme ID 087. C- 
0407. This work has been partially supported by NSF's 
grant AST-0908278. We are very grateful to Andreas 
Seifahrt for his help with LBLRTM and HITRAN, to 
Martin Kiirster for input for data analysis strategies and 
statistics, and to John Barnes for discussions about the 
search for thermal emission from hot Jupiters. We would 
further like to thank Ignas Snellen and Matteo Brogi for 
their help, and to the anonymous referee for very con- 
structive suggestions which significantly improved the 
manuscript. 



REFERENCES 



Baliunas, S. L., Henry, G. W., Donahue, R. A., Fekel, F. C, & 

Soon, W. H. 1997, ApJ, 474, L119 
Barbier-Brossat, M. & Figon, P. 2000, A&AS, 142, 217 
Barnes, J. R., Barman, T. S., Jones, H. R. A., Barber, R. J., 

Hansen, B. M. S., Prato, L., Rice, E. L., Leigh, C. J., Collier 

Cameron, A., & Pinfield, D. J. 2010, MNRAS, 401, 445 
Barnes, J. R., Barman, T. S., Jones, H. R. A., Leigh, C. J., 

Collier Cameron, A., Barber, R. J., & Pinfield, D. J. 2008, 

MNRAS, 390, 1258 
Barnes, J. R., Leigh, C. J., Jones, H. R. A., Barman, T. S., 

Pinfield, D. J., Collier Cameron, A., & Jenkins, J. S. 2007, 

MNRAS, 379, 1097 
Barrow, J. D. and Bhavsar, S. P. & Sonoda, D. H. 1982, MNRAS, 

210, 19 

Burrows, A. & Sharp, C. M. 1999, ApJ, 512, 843 

Butler, R. P., Wright, J. T., Marcy, G. W., Fischer, D. A., Vogt, 

S. S., Tinncy, C. G., Jones, H. R. A., Carter, B. D., Johnson, 

J. A., McCarthy, C, & Penny, A. J. 2006, ApJ, 646, 505 
Catala, C, Donati, J.-F., Shkolnik, E., Bohlender, D., & Alecian, 

E., 2007, MNRAS, 374, L42 
Charbonneau, D., Noyes, R. W., Korzcnnik, S. G., Nisenson, P., 

Jha, S., Vogt, S. S., & Kibrick, R. I. 1999, ApJ, 522, L145 
Clough, S. A., Iacono, M. J., & Moncet, J.-L. 1992, 

J. Geophys. Res., 97, 15761 
Collier Cameron, A., Home, K., Penny, A., & Leigh, C. 2002, 

MNRAS, 330, 187 
Cooper, C. S. & Showman, A. P. 2006, ApJ, 649, 1048 
Cubillos, P. E., Rojo, P., & Fortney, J. J. 2011, A&A, 529, A88 
Desert, J.-M., Lecavelier des Etangs, A., Hebrard, G., Sing, 

D. K., Ehrenreich, D., Ferlet, R., & Vidal-Madjar, A. 2009, 

ApJ, 699, 478 

Donati, J.-F., Moutou, C, Fares, R., Bohlender, D., Catala, C, 
Deleuil, M., Shkolnik, E., Collier Cameron, A., Jardine, M. M. 
& Walker, G. A. H. 2008, MNRAS, 385, 1179 

Gonzalez, G. 1998, A&A, 334, 221 

Grillmair, C. J., Burrows, A., Charbonneau, D., Armus, L., 

Stauffer, J., Meadows, V., van Cleve, J., von Braun, K., & 

Levine, D. 2008, Nature, 456, 767 
Gustafsson, B., Edvardsson, B., Eriksson, K., j0rgensen, U. G., 

Nordlund, A., & Plez, B. 2008, A&A, 486, 951 
Hauschildt, P. H., Baron, E., & Allard, F. 1997, ApJ, 483, 390 
Helling, C, Dehn, M., Woitke, P., & Hauschildt, P. H. 2008, ApJ, 

675, L105 

Henry, G. W., Baliunas, S. L., Donahue, R. A., Fekel, F. C, & 
Soon, W. 2000, ApJ, 531, 415 

Kaeufl, H.-U., Ballester, P., Biereichel, P., Delabre, B., 

Donaldson, R., Dorn, R., Fcdrigo, E., Finger, G., Fischer, G., 
Franza, F., Gojak, D., Huster, G., Jung, Y., Lizon, J.-L., 
Mehrgan, L., Meyer, M., Moorwood, A., Pirard, J.-F., 
Paufique, J., Pozna, E., Siebenmorgcn, R., Silbcr, A., 
Stcgmcier, J., & Wegerer, S. 2004, in Society of Photo-Optical 
Instrumentation Engineers (SPIE) Conference Series, Vol. 5492, 
Society of Photo-Optical Instrumentation Engineers (SPIE) 



Kuerster, M., Schmitt, J. H. M. M., Cutispoto, G., & Dennerl, K. 

1997, A&A, 320, 831 
Leigh, C, Collier Cameron, A., Home, K., Penny, A., & James, 

D. 2003, MNRAS, 344, 1271 
Lopez-Morales, M. & Seager, S. 2007, ApJ, 667, L191 
Press, W. H., Teukolsky, S. A., Vetterling W. T., & Flanncry 

B. P., 1992, Numerical recipes in C. The art of scientific 

computing, Cambridge University Press 
Rodler, F., Kiirster, M., & Henning, T. 2008, A&A, 485, 859 
— . 2010, A&A, 514, A23 

Rothman, L. S., Jacqucmart, D., Barbe, A., Chris Bcnncr, D., 
Birk, M., Brown, L. R., Carleer, M. R., Chackerian, C, 
Chance, K., Coudert, L. H., Dana, V., Devi, V. M., Flaud, 
J.-M., Gamache, R. R., Goldman, A., Hartmann, J.-M., Jucks, 
K. W., Maki, A. G., Mandin, J.-Y., Massie, S. T., Orphal, J., 
Perrin, A., Rinsland, C. P., Smith, M. A. H., Tennyson, J., 
Tolchenov, R. N., Toth, R. A., Vander Auwera, J., Varanasi, P., 
& Wagner, G. 2005, J. Quant. Spec. Radiat. Transf., 96, 139 

Seifahrt, A., Kaufl, H. U., Zangl, G., Bean, J. L., Richter, M. J., 
& Siebenmorgen, R. 2010, A&A, 524, All 

Sharp, C. M. & Burrows, A. 2007, ApJS, 168, 140 

Snellen, I. A. G., de Kok, R. J., de Mooij, E. J. W., & Albrecht, 
S. 2010, Nature, 465, 1049 

Swain, M. R., Vasisht, G., Tinetti, G., Bouwman, J., Chen, P., 
Yung, Y., Doming, D., & Deroo, P. 2009, ApJ, 690, L114 

Tinetti, G., Deroo, P., Swain, M. R., Griffith, C. A., Vasisht, G., 
Brown, L. R., Burke, C, & McCullough, P. 2010, ApJ, 712, 
L139 

Valenti, J. A., Butler, R. P., & Marcy, G. W. 1995, PASP, 107, 
966 

Valenti, J. A. & Fischer, D. A. 2005, ApJS, 159, 141 
van Belle, G. T. & von Braun, K. 2009, ApJ, 694, 1085 
Wiedemann, G., Deming, D., & Bjoraker, G. 2001, ApJ, 546, 1068 
Witte, S., Helling, C, Barman, T., Heidrich, N., & Hauschildt, 
P. H. 2011, A&A, 529, A44 



